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Abstract. In this paper we study the effect of a torus topology on Loop Quantum 
Cosmology. We first derive the TeichmuUer space parametrizing all possible tori using 
Thurston's theorem and construct a Hamiltonian describing the dynamics of these 
torus universes. We then compute the Ashtekar variables for a slightly simplified 
torus such that the Gauss constraint can be solved easily. We perform a canonical 
transformation so that the holomies along the edges of the torus reduce to a product 
between almost and strictly periodic functions of the new variables. The drawback of 
this transformation is that the components of the densitized triad become complicated 
functions of these variables. Nevertheless we find two ways of quantizing these 
components, which in both cases leads surprisingly to a continuous spectrum. 
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1. Introduction 

The Einstein field equations are local equations in the sense that they only describe the 
local geometry of the spacetime. For example the Robertson- Walker metric explicitly 
contains the parameter k which gives an account of the intrinsic spatial curvature. Using 
the Friedmann equations this parameter can be determined experimentally since it is 
directly related to the density parameter f2tot and the Hubble parameter h. Recent 
measurements of the energy density of the universe tend to slightly favor a positively 
curved universe [1], yet a flat curvature lies within the 1-cr range. The most direct 
conclusion is that the spatial topology of the universe is just 'E? which is the assumption 
of the ACDM model. Nevertheless in the mathematical literature it is well known that a 
flat space does not mean that its topology is necessarily R^, in fact there are 18 possible 
fiat topologies. Since the Einstein field equations are not sensitive to topology every 
possibility has to be considered as a possible candidate for the global geometry of our 
universe until it is ruled out by experiment. In order to do so we first note that the 
spectrum of the Laplace operator sensitively depends on the topology, i.e. it is discrete 
if the eigenstates are normalizable and continuous if not. In the first case the solution 
for e.g. a torus is given by plane waves with a wave vector kn taking only discrete values 
n G W while in the second case the (weak) solution to the eigenvalue equation is given 
by the (distributional) plane waves with a continuous wave vector k. For example, the 
eigenvalue problem for A on is given by (A + Efi)'^^ = 0, n G IN^, and on by 
A\l//3,/,m = (/5^ - l)\I//3,/,m, where /3elN, 0</</3-l and |m| < /. The implication of a 
solution of the form is the existence of a wave function with a maximum length 
corresponding to e.g. the length of the edges of the torus. Since the departure from 
a continuous solution is biggest for large wavelengths we have to look for large-scale 
structures of the universe in order to distinguish between cosmic topologies. The best 
way to do so is to measure the inhomogeneities of the cosmic microwave background 
(CMB), expand these in multipole moments and compare the low multipoles with the 
predictions from theory. It can be shown that in certain closed topologies a suppression 
in the power spectrum of the low multipoles is expected because of the existence of 
a largest wavelength. Since such a suppression is present in the CMB several studies 
compared the theoretical predictions for various topologies with the data. While most 
analyzed topologies can already be ruled out three of them describe the data even better 
than the infinite ACDM model, namely the torus [2], [3l H], the dodecahedron [5], [6] and 
the binary octahedron[7] (see also references therein). While the last two topologies are 
spherical the torus is the simplest model of a closed fiat topology. 

However, we know that standard cosmology cannot be the final answer as its 
predictability breaks down at the big bang. A quantization of the Friedmann equations 
a la Wheeler-DeWitt does not improve this behavior either. This situation has changed 
thanks to a new model called loop quantum cosmology (LQC) developed over the last 
few years which removes the initial singularity. LQC P[9l[l0l[IIl[12l[13l[14l[l5l[l6]is 
the approach motivated by loop quantum gravity (LQG) [171 [HI [19] to the quantization 
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of symmetric cosmological models. The usual procedure is to reduce the classical phase 
space of the full theory to a phase space with a finite number of degrees of freedom. 
The quantization of these reduced models uses the tools of the LQG and is therefore 
called LQC but it does not correspond to the cosmological sector of LQG. The results 
of LQC not only provide new insights into the quantum structure of spacetime near the 
Big-Band singularity but also remove this singularity by extending the time evolution 
to negative times. 

In sum, on the one hand we have hints from observation that our universe may 
have a closed topology, on the other hand we have a very successful loop quantization of 
various cosmologies. Thus, starting from these two motivations, we would like to study 
LQC with a torus topology. But contrary to the works on the CMB we don't want to 
restrict the analysis to a cubical torus. To do so we construct a torus using Thurston's 
theorem and find that the most general torus has six degrees of freedom which consist of 
e.g. three lengths and three angles. We will study its dynamics by numerically solving 
the Hamiltonian coupled to a scalar field. After rewritting this Hamiltonian in terms of 
Ashtekar variables we will see that the quantization of such a torus leads to a product 
between the standard Hilbert spaces of LQC and the Hilbert spaces over the circle. 
Moreover, we will find two ways to quantize the components of the triad and show that 
both (generalized) eigenfunctions are not normalizable in this Hilbert space. 

As a side remark we would like to point out that the consequences of putting 
a non-abelian gauge theory into a box with periodic boundary conditions have been 
studied in e.g. [20]. The motivation behind this idea is an attempt to explain the 
quark confinement in QCD without explicitely breaking gauge invariance. To simplify 
the analysis the sn(A^)-valued gauge field is chosen to be pure gauge, i.e. A = U^^dU 
with U G SU{N), such that the holonomy around a closed curve C only depends on 
the topological property of C. Since general relativity written in terms of Ashtekar 
variables is also a (constrained) Yang-Mills theory it may be tentalizing to use the 
methods developed for QCD in a box to LQC of a torus universe. However we will 
derive an Ashtekar connection for the homogeneous torus which is not pure gauge so 
that the holonomies along C also depend on the length of C. This may not be surprising 
in view of the fact that the Hilbert space of LQC on is spanned by almost periodic 
functions with an arbitrary length parameter fi. 

This paper is organized as follows: in Section [2] we first introduce the classical 
dynamics of a torus universe and numerically solve the Friedmann equations with a 
massless scalar field. In Section [3] we introduce the Ashtekar variables for a torus 
and also explain the complications that arise because of a closed topology. The loop 
quantization and the construction of a Hilbert space are explained in Section |4] and 
Section [5] provides a summary and directions for future works. Appendix A gives a 
short review of the fundamental domain of the 3-torus and Appendix B describes the 
dynamics of the torus in terms of Iwasawa coordinates. 
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2. Compact Homogeneous Universes and their Dynamics 

The purpose of this section is to study models in which the spatial section has a 
compact topology The compactness of a locally homogeneous space brings new degrees 
of freedom of deformations, known as Teichmiiller deformations. This leads to the 
conclusion that cosmology on a torus is simply cosmology on R'^ restricted to a cube 
may be too naive a point of view, especially since the space of solutions of a torus gets 
nine additional degrees of freedom, as already mentioned in [21]. We will introduce 
Teichmiiller spaces with an emphasis on a Thurston geometry admitting a Bianchi I 
geometry as its subgeometry [22l [23], EU [251 EE] and derive the vacuum Friedmann 
equations using the Hamiltonian formalism. 

2.1. Compact Homogeneous Spaces 

Let S be a three-dimensional, arcwise connected Riemannian manifold. 

Definition 1. A metric on a manifold S is locally homogeneous if\fp,q G E there 
exist neighborhoods U,V of p resp. q and an isometry {U,p) — > {V,q). The manifold is 
globally homogeneous if the isometry group acts transitively on the whole manifold S. 

Since S is arcwise connected we know that there is a unique universal covering 
manifold S up to diffeomorphisms with a metric given by the pullback of the metric on 
S by the covering map 

TT : S ^ S. (1) 

Singer [27] proved that the metric on S is then globally homogeneous and E is given by 
= S/F, where S is the orientation preserving isometry group of S and F its isotropy 
subgroup. 

On the other hand, we can also start from a three-dimensional, simply connected 
Riemannian manifold S which admits a compact quotient S. In order to construct 
this compact manifold consider the covering group T G S which is isomorphic to the 
fundamental group 7ri(S) of S. This implies that 

s = s/r, 

which is Hausdorff iff F is a discrete subgroup of S and a Riemannian manifold iff F 
acts freely on S. 

Definition 2. A geometry is the pair (S, S) where S a group acting transitively on S 
with compact isotropy subgroup. A geometry (S, S') is a subgeometry of (S, S) if S' is 
a subgroup of S. A geometry (S, S) is called maximal if it is not a subgeometry of any 
geometry and minimal if it does not have any subgeometry. 

We will need the following important theorem: 

Theorem 1 (Thurston [28]). Any maximal, simply connected 3-dimensional geometry 
which admits a compact quotient is equivalent to the geometry (S, 5") where S is one of 
(Euclidean), (hyperbolic). Si (3-sphere), Si x R, H'^ x R, SL{2,R), Nil or Sol. 
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If (S, S') is not a maximal geometry but is simply connected and admits a compact 
quotient as well we can find a discrete subgroub V of 5" acting freely so as to make S /V 
compact. Define (E, S) as the maximal geometry with (S, S') as its subgeometry, i.e. 
S' C S. By Thurston's Theorem (E, 5") is one of the eight Thurston geometries, which 
implies that (E, S") is a subgeometry of one of the eight Thurston geometries. 

Theorem 2. Any minimal, simply connected three-dimensional geometry is equivalent 
to (S, S), where S = R^, ^ =Bianchi I- VIII; t = S^, S =Bianchi IX; ort = Sf^xR, 
S = 5*0(3) X R, where Sp is the three-sphere and the two-sphere. 

Let Rep(E) denote the space of all discrete and faithful representations p : vri(E) — >■ 
S and the diffeomorphism : E — * E a global conformal isometry if 0*/iab = const ■ hab, 
where hab is the spatial metric of the universal covering manifold E. This allows us to 
define a relation p ~ p' in Rep(E) if there exists a conformal isometry of E connected 
to the identity with p'(a) = o p(a) o 0'. 

Definition 3. The Teichmiiller space is defined as 

Teich(E) = Rep(E)/ ~ 

with elements called Teichmiiller deformations, which are smooth and nonisometric 
deformations of the spatial metric hab ofE, leaving the universal cover {T,,hab) globally 
conformally isometric. 

The situation gets more complicated when we try to extend the previous 
construction to four-dimensional Lorentzian manifolds. The reason is that the action 
of the covering group F needs to preserve both the extrinsic curvature and the spatial 
metric of E. Thus we cannot construct a homogeneous compact manifold by the action 
of a discrete subgroup of 5* on the spatial three-section E. Instead we need the isometry 
group of the four-dimensional manifold M. Let M = R x E be a compact homogeneous 
Lorentzian manifold with metric g^j^iy and M = R x S its covering with metric g^^, 
{p,u = 0,...,4). 

Definition 4. Let (T,,hab) be a spatial section of {M^g^y). An extendible isometry is 
defined by the restriction of an isometry of (M, gab) on E which preserves E and forms 
a subgroup Esom(E) of S. 

Thus, in order to get a compact homogeneous manifold from M the covering group 
r must be a subgroup of Esom(E), i.e. 

r C Esom(E). 

The line element of M = R x E is given by 

ds'^ = -df + habit)a''a\ 

where a" are the invariant one-forms. 

Therefore the Teichmiiller parameters enlarge the parameter space by bringing new 
degrees of freedom from the deformations defined in Definition [3l In fact, the set of all 
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possible universal covers (M, cjah) carries the degrees of freedom of the local geometry 
and the covering maps F the degrees of freedom of the global geometry which are 
parameterized by the Teichmiiller parameters. 

2.2. The Torus Universe 

In this section we restrict the above analysis to the case of a flat torus and give only the 
main results. Further details can be found in [221 [23l [Ml [25l |26] . Let M = E x S be the 
universal cover of M and S the Thurston geometry (i?^, /S'0(3, R)). The isometry group 
IS0{2i) is expressed as (^(x) = Rx + a, where a is a constant vector and R e 5*0(3) in 
order that the orientation be preserved. The Killing vectors of are 

6 = dx, 6 = <9y, ^3 = d,, 

^4 = -zdy + ydz, ^5 = -xd;, + zd^,, = -yd^ + xdy. 
The line element of M is thus given by 

ds^ = -df + hahdx^'dx^ = -dt^ + a^it) ^habdx''dx\ 

where ^hab is called the fiducial metric in the LQC literature and dx"" are the invariant 
one-forms of the group IS0{3, R) \. 

* ^3 




Figure 1. The vectors oi, 02 and 03 span the torus with six Teichmiiller parameters. 
The global conformal invariance was used in order to align ai with ^1 and a2 with 
span{^i,^2}. 

The covering group F C Esom(S) = Esom(i?^) allows us to construct a torus via 
M = M/T, where M = Rx T^. The freedom of global conformal transformations allows 
us to choose the coordinate system of E such that the generators of the torus have a 
simple representation. We thus require one of the generators to be aligned with ^1 and 

I When dealing with the open case E,"^ one has to distinguish between the fiducial volume Vb of a cell 
as measured by the fiducial metric ^hab and the physical volume V as measured by the physical metric 
hab- Since we shall deal with a closed universe we have the "preferred fiducial cell" T'^ at our disposal. 
Furthermore, in the open case the spatial integrals have to be restricted to this fiducial cell whereas in 
the closed case these integrals are naturally restricted to the physical cell T'^. 
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one to lie in the ^1^2-plane. The Teichmiiller space is then generated by six Teichmiiller 
parameters in three vectors 



(2) 



where all Oj-' only depend on the coordinate time t. The configuration space C 
is therefore spanned by the six Teichmiiller parameters such that C C (see 
Appendix A ). The flat spatial metric on T"^ is then given by (a, 6=1,2, 3) 

ds^ = habdx"'dx^, hab = Oa '^(^b 
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a2 'as ' 

1\2 



a2 ^as ^ 

2\2 



(3) 



(4) 



V (sym) {as'Y + {a^Y + {as^f J 

This metric is invariant under transformations in SL^S^Ii). For example it is left 
invariant by (oi ^ — 01^,02^ —0,2', a^^ — % ^) (see Appendix A for more 



— «! , 02 ~^ ^0,2 ,a^ — ^ —o-s ) [see 
details). From Equation (j4]) we can make a Legendre transform of the Einstein-Hilbert 
action 



SE-u[g] = ^ [ *R[g], ^ = SnG 



(5) 



to obtain a Hamiltonian, where * is the Hodge star operator. After a partial integration 
of cij* (which also cancels the surface term we omitted in Equation ([5])) we find the 
following Lagrangian: 

1 1 



X 



+ (02 



Ak ai ^02 ^03 ^ 

((«3 r + (as r) («2 'a, ' - a, 'd2 ') ' + (a, 'Yias 'f{d2 J 

2\2 



2ai 'as 'di 'd^ ' 



+ ai ' (^ai ' ((ds + (d3 - 4a3 ^di 'd-s 
2a2 ^|a2 ^03 ^di ' (03 ^di ' — ai 'ds ') 



+ ai 



ai 'as ^d2 ^d3 ' 



as 'as ^di ^d2 ' 



+ d2 ^ (ai ^03 ^d3 ^ + 2a3 ^(03 ^di ' + ai 'd^ ^)) 



We introduce the momenta 
'^-■=da? 



(6) 



conjugate to the configuration variables such that the phase space V = T*C C R 
is the cotangent bundle over C with 

{aa '',p'd} = 5l5i {aa ^ a, ^} = 0, p'^ = 0, 



12 



(7) 
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where the Poisson brackets are defined as 

3 



dtta^ dp"" b daa ^ dp"" b 

a, 0=1 

for any smooth functions on the phase space. We insert (ia * = Oa ''{p'^ d) into the Legendre 
transform of Equation ([5]) and get the Hamiltonian 

(ai y 1)^ + {a2 V 2? + (a3 V 3)^ + [a, V 0^ + 4(a2 V 1)^ 
+ y 1)^ + 4(a3 y 1)^ + 4(a3 y 1)^ + (a3 y 2)^ 

+ 4(03 y 2)' -203^03 y2p'3 

+ 2ai 1 (^2 1 - 02 y 2 + 03 1 - as y 2 - 03 y 3) (S) 

- 203 1 (03 2 + 03 y 3) 

- 2a2 1 (^2 2 - ^3 + as 2 + as s) 

+ 2a2 ^jaa ^ (4^^ ^ + p2 ^^s 2) _ ^2 ^ (ag y ^ + 03 y 3) } 

The Hamiltonian constraint Tig ~ reduces the dynamical degrees of freedom from 
dim P = 12 to dim V = 10, which agrees with [21]. To compare this Hamiltonian with 
the usual Bianchi type I models we set all offdiagonal elements to zero and aj* = a^, 
p\ = p^ (no summation) , and get 

^ 4 \ 0203 aia^i aia2 a-^ ai a2 / 

which agrees with the result given in [31] up to a factor 2 in the definition of the action. 
To get the isotropic case§ we further set ai = a, = p/3 and find that the Hamiltonian 
([8]) reduces to the usual first Friedmann equation 

rio 



12a 

and the Hamiltonian equation p'^ j = —dTig/dai^ to the usual second Friedmann 
equation 

dHg np^ 



P- ' 



da 12a2 ' 

The second Hamiltonian equation is given by d = dHg/ dp = —up/ (6a) and allows us to 
recast the first Friedmann equation into the usual form 

1-Lg = —'iaa^ I K. 

Furthermore, notice that all ai^ and j, i 7^ j, have to vanish in order for the torus to 
remain aligned with the Killing fields ^7. 

§ At this point care has to be taken because there is no homogeneous and isotropic vacuum solution 
to the Einstein equation (see Section [3.4. ip . Only a nonvanishing energy-momentum tensor allows for 
the isotropic limit of Tig, which corresponds then to the usual Friedmann solutions. 
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Figure 2. Left panel: Solutions corresponding to the Hamiltonian pTT)) with the 
initial conditions ai'(l) = 1, = —1 (no summation), ai^(l) = (i 7^ j), 

p^i{l) = ^^2(1) = 0, 1(1) = 0.2, 0(1) = 10-3, 7r(l) = 1.2. The diagonal momenta 
i are chosen to be negative such that all sides of the torus expand. The solid black 
line is ai ^, the dashed one 02 ^, the dotted one 03 ^ and the gray one the ofF-diagonal 
02 ^ . The time t parametrizes the coordinate time in natural units (c = K = ?i=l). 
Right panel: Solution corresponding to the Hamiltonian (jlOp at two different times. 
The initial condition is a cubic universe with a;* = ao, p^ i = po, = (i 7^ j), 
P'i 7^ (i 7^ j). For both panels the mass and the potential of the scalar field have 
been set to zero. 



We add a matter term consisting of a homogeneous massive scalar field|| to Equation 
(IH]) to obtain the Hamiltonian 

H = Hg + n^ = n, + ^TT^ + ^m|02 + v^y(0), (10) 

where tt is the momentum of the scalar field, its mass, h = (ai ^)^(a2 ^)^(a3 ^)^ the 
determinant of the spatial metric (j4]) and V{(f)) the potential which we set to zero in 
the sequel. From this equation we calculate the Friedmann equations and compute 
the shape of the universe for a special choice of initial conditions, which is shown in 
Figure O All classical solutions have the limit limt^o = and grow with a oc 
for a massless scalar field with zero potential. Furthermore, note the convergence of ai ^ 



and a2 ^, which is explained in [Appendix B 



3. Symmetry Reduction and Classical Phase Space for Ashtekar Variables 

In this section we shall repeat the complete analysis introduced in [291 [9l [lO] in order to 
see the role of a compact topology on a connection. Our strategy is to find an invariant 
connection on the covering space M and then restrict it to the compact space M by 

II Notice that since every scalar field lives in the trivial representation of the rotation group it is not 
possible to construct a scalar field which is homogeneous but anisotropic. 
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means of the covering map ([T]). In the following section, when referring to the covering 
space, we shall use a tilde. 

3.1. Invariant Connections 

Let P{M, SU{2), it) be a principal fiber bundle over M with structure group SU{2) and 
projection tt : P ^ M. We require that there be a symmetry group S C Aut(P) of 
bundle automorphisms which acts transitively. Furthermore, for Bianchi I models S 
does not have a non-trivial isotropy subgroup F so that the base manifold is isomorphic 
to the symmetry group 5', i.e. M/S = {xq} is represented by a single point that 
can be chosen arbitrarily in M. Since the isotropy group F is trivial the coset space 
S/F = S is reductive with a decomposition of the Lie algebra of S according to 
CS = CF ® CF± = CF^ together with the trivial condition Ad^^F^ C CF^. This 
allows us to use the general framework described in [9], [TOl [29] . 

Since the isotropy group plays an important role in classifying symmetric bundles 
and invariant connections we describe the general case of a general isotropy group F. 
Fixing a point x G M, the action of F yields a map F : tt~^{x) tt~^{x) of the 
fiber over x. To each point p G tt~^{x) in the fiber we assign a group homomorphism 
Xp : F ^ G defined by f{p) =: p ■ Xp{f), V/ G F. As this homomorphism transforms 
by conjugation Xp.g = Ad^-i o Ap only the conjugacy class [A] of a given homomorphism 
matters. In fact, it can be shown [29] that an S'-symmetric principal bundle P{M, G, vr) 
with isotropy subgroup F C 5 is uniquely characterized by a conjugacy class [A] of 
homomorphisms \ : F ^ G together with a reduced bundle Q{M /S, Zo{X{F)),'71q), 
where Zg{X{F)) is the centralizer of X{F) in G. In our case, since F = {1} all 
homomorphisms X : F G = SU{2) are given by 1 i— l^. 

After having classified the S'-symmetric fiber bundle P we seek a [A] -invariant 
connection on P. We use the following general result [30] : 

Theorem 3 (Generalized Wang theorem). Let P be an S -symmetric principal bundle 
classified by ([A], Q) and let u be a connection in P which is invariant under the action 
of S. Then u is classified by a connection uq in Q and a scalar field (usually called the 
Higgs field) (p : Q x CF± —>■ CG obeying the condition 

0(Ad/(X)) = AdA(/)0(X), V/ G F, X G CF^. (11) 

The connection uj can be reconstructed from its classsifying structure as follows. 
According to the decomposition M = M/S x S/F we have u = uiq + u^s/f with 
^s/F = ° '•*6'mc5 where l : S/F ^ 5* is a local embedding and 9mc is the Maurer- 
Cartan form on S. The structure group G acts on by conjugation, whereas the 
solution space of Equation ( iTTl) is only invariant with respect to the reduced structure 
group Zg{X{F)). This fact leads to a partial gauge fixing since the connection form 
ojg/p is a ZG(A(F))-connection which explicitly depends on A. We then break down the 
structure group from G to Zg{X{F)) by fixing a A G [A]. 
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In our case, the embedding l : S ^ S i?, the identity and the base manifold 
M/S = {xq} of the orbit bundle is represented by a single point so that the invariant 
connection is given by 

A = (j)o Que- 

The three generators of CS are given by Tj, 1 < / < 3, with the relation [Tj,Tj] = 
for Bianchi I models. The Maurer-Cartan form is given by 6mc = Cj^Ti where uj^ are 
the left invariant one-forms on 5*. The condition ffTTj) is empty so that the Higgs field 
is given by : CS CG, Tj i— s> 0(T/) =: (pi^'Ti, where the matrices tj = — icrj/2, 
1 < j < 3, generate CG, where aj are the standard Pauli matrices^. In summary the 
invariant connection is given by 

A = (piWidu^. (12) 

In order to restrict this invariant connection we define the invariant connection A on 
with the pullback given by the covering map ([I]). The generators of the Teichmiiller 
space (see Equation ([2])) allow us to write A as: 



I A. 1 
V 








02 





(13) 



3.1.1. Simplified Model In the sequel we shall concentrate on a simpler model for 
which we can also easily satisfy the Gauss constraint. We choose a torus generated by 
the vectors ai = (ai^,0,0)"^, a2 = (0,02^,02^)"^ and = (0,03^,03^)-^ (see Figure [3]) 
such that 



( ^1' 


V 



\ 



and 







/ ai^ \ 

a2 ^ 03 ^ 

y 02^ 03^ / 

\ 



02^02^+02^03^ 03^02^+03^03- 



V 



O2 ^03 2 
02 V2 ^ + O2 ^03 ^ 



O3 ^2 ^ + O3 303 ^ J 



The vectors Xj dual to are given by 







/ \ 



03' 



V / 

where we defined f) = 02 ^03 ^ 
% We use the convention TiTj = ^Cij ^Tk — jSijt 2x2 



X. 



( \ 

-03^ 

V J 



02^03 ^ 



(14) 
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3.2. Classical Phase Space for Ashtekar Variables 

The phase space of full general relativity in the Ashtekar representation is spanned by 
the SU (2)-connection A\ = V\ + 7/^* and the densitized triad = \ det e|e^, where 
is the spin connection, /T* the extrinsic curvature, the triad and 7 > the Immirzi 
parameter [TTl [TBI [19] . The symplectic stucture of full general relativity is given by the 
Poisson bracket 



{Al{y),E'^ix)} = ^6t6}6{x,y). 
The connection between the metric and the densitized triad is given by 



hh 



ab 



d'^EfE], 



(15) 
(16) 



where h is the inverse of the metric hat- 

We can now use the results obtained in Section [2] to construct the phase space V in 
this representation. In the preceding subsection we have already found the configuration 



space is spanned by (see Equation (IT^ 
dual to the connection is given by 

/ 

/ ail 

{E^) = ^p^X^j = ^/h 

V 



On the other hand, the densitized triad 



13 P 2—a-i P 2 03 P 3—03 P 3 



. 2i;3 



. 3^;2 



0.2 "P" 2-0.2 "P" 2 0.2 ^P^ 3 -"2 '^P' 3 



, 3^;2, 



(17) 



where 



I P^i \ 



together with uj^-^i 



\ 








P^2 P% 

P^3 / 



6 J and h 



l\2l 



3\2 



— aiiu. lu — {ai"y{a2"a^" — a2"a^")" is the determinant of 
the spatial metric constructed from the vectors and j the momentum dual to 0/ * 
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satisfying the Poisson bracket 

{h\f ,} = ^^i8] (18) 

with the volume Vq = d^x\/h of as measured by the metric h. For later purpose 
we define new variables 

0/ = L,0/, p', = ^p'„ (19) 

{<Pj\p',} = kjS^5}, (20) 



such that 



where 

Li = a,\ L, = v/M^T(VF, L, = ^{a^Y + ia^^y. 
Thus we conclude that 

Proposition 1. The classical configuration space As = is spanned by the five 
configuration variables 0/*. The phase space V = M.^^ is spanned by 0/* and the five 
momenta p"^ j satisfying the Poisson bracket (d^. 

Furthermore, note that the determinant of the densitized triad is given by 

det = tp' i(/ 2 - 2P^ 3), (21) 

where we defined 



Vr 







The relation between the new variables {ipi^iP"^ j) and the 'scale factors' a^'' and 
their respective momenta p"" h can be found by using Equation ( fT6l) and the Poisson 
brackets ([7]) and (!20|) . A closed form could only be found for p^ 1 and is given by 

bi I I 2 3 3 2 1 
l| = 1^2 03 - 02 03 |. 

3.3. Constraints in Ashtekar Variables on the Torus 

In the canonical variables (IT^ the Legendre transform of the Einstein-Hilbert action 
([5]) results in a fully constrained system [T71 [HI [19] 

S= — [ dt [ d^x (2A^aEt - WGj + N^Ha + NH]) , (22) 
2/^ J-R Jt^ ^ ' 

where Gj is the Gauss constraint, Ha the diffeomorphism (or vector) constraint, 

the Hamiltonian and A-', A^", A^ are Lagrange multipliers. The Hamiltonian constraint 

simplifies to 

Cgrav = -- / rf^xA^e,,,F^,^== (23) 

due to spatial flatness, where the curvature of the Ashtekar connection is given by 

^ab = daAl — dbA''^ + e'ji^A{A'^ = e''ji^A{A'^. 
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Homogeneity further requires that 7^ N{x). Inserting Equation fll4p and Equation 
f|T7|) into Equation (123|) we get 

1 1 



grav 



+ 



'<P^-<i>2'<P,'){p\p' 



P 2P 3) 



(24) 



where we defined N = L1L2L3/V0 in order to simphfy the Hamiltonian. Using the 
Hamiltonian (l24l) we can compute the time evolution of the basic variables 0i ^ and j 
(see Figure H]). Setting all off-diagonal terms to zero we see that Equation ( l24l) matches 
with Eq. (3.20) in [31]. If we further set 0(j) * = c and p^*-* , = p we get 



grav 



which is exactly the same result as the homogeneous and isotropic case [S]. 




1.02 1.04 1.06 



Figure 4. Solutions corresponding to the Hamiltonian (|24p coupled to a massless 
scalar field with vanishing potential. Left panel: the black thick solid shows the 
evolution of 0i ^, 4)2 '^ is the black dashed line, 03^ the dotted line, 02'^ the gray 
dashdotted one and 03 ^ the solid gray one. Right panel: the black thick solid shows 
the evolution of 1, 2 is the black dashed line, p^ 3 the dotted line, p^ 3 the gray 
dashdotted one and p^ 2 the solid gray one. In both cases the initial conditions are 
01 1 = 1.0, 02 2 = 0.2, 03 3 = 0.4, 02 3 = 0.6, 03' = 0.7, pS = 1.0, 2 ^ 0.3, 
p^ 3 — 0.5, p^30.5, p^ 2 — 1.4, = 0.01 and = 8.1. The time t parametrizes the 
coordinate time in natural units {c ~ k = h = 1). 



3.4- Diffeomeorphism and Gauss Constraints 

The Gauss constraint stems from the fact that we chose the densitized triads as 
the momenta conjugated to the connections A^^. In fact, the spacial metric can be 
directly obtained from the densitized triads through Equation (fT6l) and is invariant 
under rotations given by Ef k-s- OlEj. In order that the theory be invariant under such 
rotations the Gauss constraint 

Gi = daE^ + e^.uKEl ^ (25) 
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must be satisfied. The diffeomorphism constraint modulo Gauss constraint originates 
from tlie requirement of independence from any spatial coordinate system or background 
and is given by 

= Fi,E\ ^ 0. (26) 

However, as mentioned in [21] we have to be careful when dealing with these 
constraints in the case where the topology is closed. We thus divide this subsection 
into two parts, starting with the general case of open models. 

3. 4-1- Open Models Due to spatial homogeneity of Bianchi type I models the basic 
variables can be diagonalized to [HI [Hj 

where u' is the left-invariant 1-form, X' the densitized left-invariant vector field dual 
to uj' and A' G 5*0(3)"'". This choice of variables automatically satisfies the vector 
and Gauss constraints, thus reducing the analysis of Equation fl22|) to the Hamiltonian 
constraint ( l23l) . The homogeneous, anisotropic vacuum solution to the Einstein field 
equations is called the Kasner solution and is given by the following metric: 

ds^ = -dr^ + r^'^'dxl + r^^^dxl + r^'^'^dxl 

where the two constraints ^ R, = = 1 have to be fulfilled. These imply 

that not all Kasner exponents can be equal, i.e. isotropic expansion or contraction 
of space is not allowed. By contrast the RW metric is able to expand or contract 
isotropically because of the presence of matter. At the end, from the twelve-dimensional 
phase space only two degrees of freedom remains. 

An infinitesimal diffeomorphism generated by a vector field V induces the following 
action on the left-invariant 1-form a;': 

u;a^u;'a + eCyi^a, (27) 

where Cy is the Lie derivative along V. Such transformations leave the metric 
homogeneous provided the vector fields satisfy 

= -uy)K'' (28) 

for some constants /j and functions given by i^KjU^ = ^] [21]. The last equation for 
relies on the fact that the 3-surface is topologically R'^ and the Killing vectors Ki 
commute. As we shall see below this will not be the case in the closed models. 

In the case of rotational symmetry the diffeomorphism constraint is once again 
satisfied by the choice of variables whereas the Gauss constraint is not. However, in 
such a case the triad components can be rotated until the Gauss constraint is also 
satisfied. Further details can be found in |32j . 

"*" In order to avoid confusion with the rest of this work we tag every variable with a ""' when deahng 
with the open case. 
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3.4-2. The Torus as a Closed Model As we have seen in Section |2] it is not possible 
to align the Killing fields with the left-invariant vectors, whence the metric takes the 
non-diagonal form (j4]) and the Ashtekar connection the form ( fT3l) . In the previous 
subsection we saw that a diffeomorphism preserves homogeneity provided it satisfies the 
condition ( l28l) . In the closed model the analysis goes through as well and we find that 
Vi has to satisfy the same condition fl28|) . However, since such fields lack the required 
periodicity in we are led to the conclusion that there are no globally defined, non- 
trivial homogeneity preserving diffeomorphisms (HPDs) and there is no analog of (1271) . 
Thus, instead of one degree of freedom we get additional degrees of freedom. 
The Gauss constraint for a Bianchi type I model is given by 

Gi = €ijk<PiV k. (29) 

With our choice of variables two Gauss constraints are automatically satisfied, namely 
^2 = ^3 = 0. However, we can still perform a global SU{2) transformation along ti 
which is implemented in the nonvanishing Gauss constraint 

Gl = 02 V 3 + hV3- 02 V 2 - 03 'p' 2 ~ (30) 

generating simultaneous rotations of the pairs (02^,02'^), (p^2,P^3) resp. (03^,03^), 
{p^2jP^ 3)- Thus the norms of these vectors and the scalar products between them are 
gauge invariant. The Gauss constraint allows us to get rid of e.g. the pair (03^,p'^2) 
by fixing the gauge in the following way: we rotate the connection components such 
that 03^ = 0. Because the length ||03|| = v^^03^P~-r"(03^)^ is preserved we know that 
03^ 7^ 0. The Gauss constraint then implies that 2 = (02 3 — 02^2)703^- This 
gauge fixing reduces the degrees of freedom by two units. 

The diffeomorphism constraint is given by Equation (12^ and since F^f^ = e* jkA^^A\ 
{daA\ = thanks to homogeneity) we find that 

= ,,AiA',E^ ^ Ai^G,. (31) 

The gauge fixing we just performed ensures that the diffeomorphism constraint also 
vanishes. 



3.5. Canonical Transformation 

In this subsection we introduce a set of new variables which will greatly simplify the 
analysis of the kinematical Hilbert space. We first perform a canonical transformation 
on the unreduced phase space: 

Qi = 0i\ P'=p\, 

n _,/7X^\^— TX^ p2_ pS02'+P%02' 

"^^-^^^^^ "V(02T + (02T 

Q3 = Vih'y + ih')', P' = P '^'^^'^P'''^'' . (32) 
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2^2 



+ 



3^2 



02 



-P^ 303 ^ + 203 ^ 



such that the variables are mutually conjugate: 



We choose the convention that the diagonal limit can be recovered by setting 6i = 62 = 0. 
The inverse of this canonical transformation will be important in the sequel and is given 
by: 



^2 



P^2 



Q2 COS 

p2 cos(^i 
Q3 sin 
P^ sin I 



1^ 



Pdi sin( 



'2), 
'2) + 



Q2 

Pg, C0S(6'2) 



Q3 



Q2sin(6'i), 

Pa, COS 



C0S(6'2), 

Pe^ sin (6*2 



i^ + p2sin(^i), (33) 



Q3 



P'^ cosf 



It is important to note that Q21 Q3 € R+ and 6*1, ^2 £ [k7T, {k + l)^] where we restrict the 
values of k to be either A; = if sgn(02 > or /c = 1 if sgn(02 ^) < 0. If sgn(02 = 
then we have the case /c = if sgn(02 ^) > or A; = 1 if sgn(02 ^) < 0. The function 
arcicos(x) is related to the principal value via arcicos(x) = 27r— arccos(x). With this 
convention we can recover Equation 0321) unambiguiously from Equation (1331) . 
The Hamiltonian constraint fl2^ is given in terms of the new variables by 

(2/€72)-i 



a 



grav 



X 



Pl[cos(ei+e2)(Pfl^Pfl^-P2p3Q2Q3) + (i"^i"92Q2+P3PfllQ3)sm(ei+e2) 

Q2Q3 



xhp^Qi cos(2^2)P^?2 + P'Q2(sin(2^i) 



P3g3(sin(2^2 



P^Q2 Pe, sin(2(^i + 62)) -2 cos^Oi + 92)P^Qs 



(34) 



+Pe, 



2 cos2(^i + e2)Pe, + 2 cos(2^i)pigi + P'^Pe, sin(2(^i + ^2)) 



Using this Hamiltonian we can compute the time evolution of the basic variables Qi, 
6a, P* and Pg^ (see Figure [5]). We choose the initial conditions so that they correspond 
to the values of the old variables (see caption of Figure Hj). By doing so we are able to 
check whether the solutions to Equation ( 12^ are equivalent to the solutions to Equation 
flM|) by performing the canonical transformation fl32|) . The different solutions do indeed 
match up to a very good accuracy. 

The only nontrivial Gauss constraint (130|) is then given by 

Gi = Pg^ — Pg^, (35) 

which vanishes only when Pg^ = Pg^. We are free to fix the gauge by setting 62 = 0. The 
same result can be obtained from the gauge fixing performed in Section 13.4.21 so that 



g3 = 03', ^^' = ^'3, 02 = and Pe,=Pg,. 
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I.OO 1.02 1.04 1.06 1.08 1.10 1. 1.02 1.04 1.06 1.08 1.1 

t t 



Figure 5. Solutions corresponding to the Hamiltonian (j34p coupled to a massless 
scalar field with vanishing potential. Left panel: the black thick solid shows the 
evolution of Qi, Q2 is the black dashed line, Qs the dotted line, 9i the gray dashdotted 
one and 02 the solid gray one. Right panel: the black thick solid shows the evolution 
of P^, is the black dashed line, the dotted line, Pe^ the gray dashdotted one 
and the solid gray one. In both cases the initial conditions arc Qi = 1, Q2 = 0.63, 
= 0.81, = 1.25, 92 = 1.05, P^ ^ 1, P^ = 0.57, P^ = 1.46, Pg, = -0.08, 
Pe^ — 0.21, (f) ~ 0.01 and = 8.1. The time t parametrizes the coordinate time in 
natural units {c = k = h = 1). 

The symplectic structure of the reduced 8-dimensional phase space is given by 
n = —(dQi A dP^ + dQ, A dP'^ + dQ3 A dP^ + dOi A dPe,). 

4. Kinematical Hilbert Space 
4.I. Holonomies 

In the last section we parametrized the classical phase space and gave the Hamiltonian 
in terms of the Ashtekar variables. To quantize the theory we have to select a set of 
elementary observables which have unambiguous operator analogs. In order to do so we 
have first to find elementary variables of the 5-dimensional configuration space. 

According to [TU] the configuration space on the covering space E is given by Higgs 
fields in a single point xq which is the only point in the reduced manifold S/S*. In 
quantum theory these fields are represented as point holonomies associated to the point 
Xq On S we can take the three edges ^1, .^2 and .^3 in order to regularize the point 
holonomies. However, we would like to apply this construction to a closed manifold. 
First we note that H = S /T such that two elements g,g' ^ S are equivalent if there 
is an element 7 G F such that g' = g + We can thus restrict the regularization of 
the point holomonies to the three edges Xi, X2 and X3 meeting at xq without losing 
information. Our elementary configuration variables are then the holonomies along 
straight lines 77 : [0, 1] ^ E defined by the connection A oc {(j)/) [TOl [TTl [T2l [T3l [T4] . 
Now the holonomies along Xi, X2 resp. X3 are given by 

h^^^^ = exp(Ai0i ^Ti) = cos(Ai0i ^) + 2ri sin(Ai0i ^) 
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h'i''^ = exp(A202 V2 + A202 V3) (36) 

/l("^)= exp(A303V2 + A303V3), 

where A/ G (— cxd, 00) and XjLj is the length of the edge / with respect to the spatial 
metric hat- The auxilary Hilbert space is then generated by spin networks associated 
with graphs consisting of the three edges 7/ meeting at the vertex xq- 

In open Bianchi type I models the gauge invariant information of the connection 
can be separated from the gauge degrees of freedom via the relation (p] = C(/)A} with 
A G 5*0(3) so that the holonomies become simple trigonometric functions. In our case 
the situation is more complicated because the holonomies h2 and cannot be reduced 
to such functions since 

= cos (Ac, II 4 II 72) + sin (x^\\$^\\/2) (no summation), 

where a = 2, 3 and 




The problem is that this expression cannot be used in this form since there is no 
well defined operator 0/ * on the kinematical Hilbert space. Using the canonical 
transformation fl32l) we can re-express the holonomies such that 

h^^'^ = cos(Aigi/2) + 2ri sin(Aigi/2), 

hi^'^ = cos(A2Q2/2) + 2(r2 cos^i + r-, sin ^i) sin(A2g2/2), (37) 

= cos(A3g3/2) + 2(r2sin^2 + r3Cos^2) sin(A3g3/2). 

Since A/ G R matrix elements of the exponentials of Qi, Q2 and g3 form a C*-algebra 
of almost periodic functions. On the other hand the variables 61^2 are periodic angles 
such that only strictly periodic functions exp(ika6a) ^ ^(1) with /c^ G Z are allowed. 
Thus, any function generated by this set can be written as 

Al,A2,A3,fci,fc2 

X exp ^iXiQi + ^i\2Q2 + ^ihQs + ihdi + ik26^ (38) 

with coefficients ^Ai,A2,A3,fei,fc2 ^ C, generating the C*-algebra As- Note that this function 
is almost periodic in Qi,Q2 and Q3 and strictly periodic in 61 and 62- The spectrum 
of the algebra of the almost periodic functions is called the Bohr compactification 
Rb '■= A(Cyl5) of the real line and can be seen as the space of generalized connections 
[T3| . Thus the functions (l38l) provide us a complete set of continuous functions on 
X R^ X Rb X X S^. Moreover the Gel'fand theory guarantees that the space R^ 
is compact and Hausdorff [35] with a unique normalized Haar measure dfi{c) such that 

1 

f{c)d^i{c):= \mv^— f{c)dc. 
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A Cauchy completion leads to a Hilbert space Ti^ defined by the tensor product 
= nf ® Hi? with the Hilbert spaces Kb = L'^{RB,dfi{c)) and Hsi = L'^{S\d(l)) 
of square integrable functions on and the circle respectively, where d(j) is the 
Haar measure for S^. An orthonormal basis for Ti^ is given by the almost periodic 
functions {Qi\fii) = exp(z/i/Q//2) (no summation) with yU/ e R with = S^^y^. 

Analogously a basis for Ti^i is given by the strictly periodic functions {Oa\ka) = 
exp{ika9a) with {K\k'J = 4,.,fc^,- 

We choose a representation where the configuration variables, now promoted to 
operators, act by multiplication via: 

{gi92){Q,9) = gi{Q,9)g2{Q,9). 

The momentum operators act by derivation in the following way: 

The eigenstates of all momentum operators are given by 

l/i, k) := |/ii,/i2,/i3, ki, k2) 

:= l/ii) (g) l/ia) ® I/U3) ® \ki) ® \k2) 

with 

P'lfl, k) = lllyfillfl, k), k) = ^k^lfl, k). (40) 

The simple form of the momentum operators ( l39l) may suggest that the Hilbert 
space of LQC on a torus is simply expanded from L^(R|) to L^(R|) x L^([/(l)^). 
However the situation is far more complicated because the important variables for 
the Gauss and Hamiltonian constraints are not the new momenta and Pq but 
the components i of the triad. In terms of the new canonical variables they are 
complicated functions of both the configuration and momentum variables, as can be 
seen from Equation (!33|) . These expressions cannot be quantized directly since the 
operators Q2,3 fail to be well defined on the Hilbert space. The solution is to consider 
the momentum operators of the full theory given by a sum of left and right invariant 
vector fields. In [TT] the same strategy was used to show that the triad components 
i act by derivation. In our case the situation is more complicated since the triad 
components contain both configuration and momentum variables. The triad operators 
act on functions in Ti^ and are given by 

P^ = -^(^f^(M+xl'^(M), (41) 

where xf\hi) and xl^\hi) are the right and left invariant vector fields acting on the 
copy of SU (2) associated with the edge e/ of length 1 and are given by 



xf)(M=tr 



d 

dhi 



xP(M=tr 
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Applying the operators 2 and 3 on the function tr(/i2) we get 

p2 2 tr(/i2) = 2^^ 2 cos(A2(52/2) = z7^piA2 sin(A2(52/2) cos(^i), 

p2 3 tr(/i2) = 3 cos(A2g2/2) = z7^^iA2 sin(A2Q2/2) sin(^i). (42) 

We see that the usual expressions for an open topology can be recovered by simply 
setting 61 = 0. Applying these operators once again we get the expressions: 

which means that cos(A2Q2/2) is an eigenfunction of both {p^ 2Y and (p^s)^ with 
eigenvalue 7^/pj/2A|. On the other hand we have 

P^ 2P^ 3 tr(/i2) = p^ 2 tr(/i2) = 0. 
4-2. Quantization: 1. Possibility 

As previously mentioned we cannot directly quantize the expressions f l33|) because 
Qi does not exist as multiplication operator on Ti^ . In a loop quantization only 
holonomies of the connections are represented as well-defined operators on Ti^ . Thus we 
replace every configuration variable Qi in Equation (!33|) by sin(5/(5//2)/5/ [36], where 
5i G 1EI\{0} plays the role of a regulator, and compare it with the results just obtained in 
terms of left and right invariant vector fields. For later purpose we order the operators 
in a symmetrical way get the following operators acting on functions of Ti^: 



^2' 



^3 



2 Sm 52Q2 n -2 n r>2 '^2Vsin^i ^ r— - 

= 7 cos6'i, p 2 = cos6'iP - ^— — — Vsmfc'i, 

02 sm(d2Q2) 

sin(52Q2) . a .2 • /) £>2 , 52Vcos^i ^ I — 

sma*!, p 3 = sm6'iP + — — — — PeiVcos6'i, 

02 sm(d2Q2) 

2 sin(53Q3) .3 • /) p3 , 53Vcos^2 A / ^ 

= 1 sm6'2, p 2 = sm6'2P + ^— — — PeaVCOsa'a, (43) 

03 miyp-iQ-i) 



I3' 



sin(^3Q3) . .3 _^^,z,p3 ^3Vsing2 A ^/— 7- 

COSy2, P 3 — COSy2-r — -Tea V 8111^2. 

t'3 sm(^3g3) 



Applying e.g. the operator 2 on cos(A2(52/2) with the definitions PUI) we see that we 
obtain the same result as Equation ( l42i) for 5 = 1. This is not surprising in view of the 
fact that we defined the operator p^ j in Equation (14T|) with holonomies along edges e/ 
of length 1. 

This substitution is problematic since the configuration variables (^2,3 are by 
definition positive (see Equation fl32|) ). Therefore, for (^2,3 sin(52,3Q2,3)/(^2,3 to be 
valid we restrict the analysis to the domain < (^2,3 < tt. In the diagonal case the 
situation is less problematic because the configuration variable c is arbitrary such that 
sin(5c) is also allowed to be negative. 

Classically, since the change of variables (l33l) is a canonical transformation the 
symplectic structure is conserved, i.e. the Poisson bracket between 2 and p'^ 3 vanishes: 

{p^2,p\}qp = 
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A quantization of the above expression is obtained with the substitution {, } — "^f, ]^ 
such that the commutator between 2 and p'^ 3 should also vanish. However, the 
consequence of the substitution of 1/Qi by 61/ sm{6iQi) is that commutator between 
these two variables doesn't vanish anymore: 

[-^2 -2 1 r^n Q\ 2,4 . C0S(52<52) - ^2 5/ 

sm^(d)2Q2) 061 

Formally we can recover the classical limit by taking the limit 

lim[p\,p\]fiQ2,e^) = 0, 
02— >o 

which however fails to exist on Ti^ . 

The operators p^ i are partial differential operators with periodic coefficients in 
both 6 and Q. In spherically symmetric quantum geometry a similar situation arises 
when considering the quantization of a nondiagonal triad component [36]. However the 
expression of this component reduces to a Hamiltonian whose eigenvalues are discrete. 
In our case the situation is more complicated. 

4-2.1. Quantization of p^ 2 In order to find eigenfunctions of the triad operators let us 
consider an operator of the form 



■= -zcosg— + 2 . — Vsmg. 

A substitution ^ = 6Q shows that As = 6A1 = 6A so that it is sufficient to determine 
the spectrum for 6=1. This operator is symmetric on Ha '■= L^{^b, dfis) ® -^^(^(1))- 

{f,Ag) = {Af,g), yf,geV{A), 

where T>{A) C Ha is the domain of A. The eigenfunctions of A are obtained by solving 

Af^{^,e) = xm,e),i.e. 

. sine z cos ^ 

sm t, oO 2 sm 4 

where we constrain ^ to be in the interval [0, vr] in order to avoid negative values of sin,^. 
We look for a solution of the form w = w{^, 6) [37] satisfying 

dw sin 6 dw ( i cos 9\ dw 

—I cos 0-— + I— TTTTT = 1 ~ 7;~ 7 /A" 



sin do \ 2 sin ^ / df\ 
such that the characteristic functions are given by 

smut) V 2sm^(t)y 

where the dot is the time derivative. Combining the first two equations gives after 
integration 

sinetan| = Ci, (47) 
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meaning that every C^-function Oi(sin0tan(^/2)) solves the left-hand side of Equation 
l5|) . In order to solve Equation fH5l) for A 7^ we first note that 



cose(t) = cos(arcsin(Ci cot(^/2))) = ^Jl-Cf cot2(^/2) = i^. (48) 
An integration of this equation gives the result 

^ ^ , x/261og(v/2acos(^/2) + &) 
'aVl-Ci^cot^(^/2)|sin(e/2)|' 
where 



a=^Jl + Cl and b = ^J -1 + Cf + cos^{l + Cf) 
The last characteristic equation in P6l) can be written as 



V 2sin{ 

such that 



9.^ \ COS 9 2 sin ^ ^ 
Equation (HHl) can be inserted into the last equation such that after an integration we 
get the result 

log fx = Xt + log (v/tan(e/2)) + C, 



where t is given by Equation (1491) and C is an integration constant. The final solution 
to the PDE (jinD is thus given by 

/A(e,^)=Mv/tan(e/2)x 



(^V2cos(e/2)ai + /5i) fii(sin^tan(e/2)), (50) 



ai{^,e) = a/1 + sin2 0tan2(^/2) and 



where 



p^{^,9) = y'-l +cos^+ (1 + cosOtan2(^/2)sin2^. 

The C^-function f2i(sin^tan(,^/2)) can be determined by e.g. boundary conditions. For 
simplicity we set f2i(sin6'tan(^/2)) = 1 subsequently. As a cross-check we see that the 
first line of Equation (150|) solves 



^9vtan(£/2) i cos9 , ,^ , , 

-zcosg ^ ^ + -— -v/tan e/2 = 

(74 2 sm 4 

and the second one the eigenvalue problem of the operator A. The scalar product on 
a. A is given by 

(A,/v)= hm ^\ di\ dOhh. (51) 

The integral of | -^tan(4'/2)p over one period is not finite and since the second line of 
Equation (1501) never vanishes the eigenfunctions fx are not normalizable in 7i^. We 
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could choose the function Vti oc (tan(^/2))~^/^ but we would automatically get the 
factor (sin 6')"^/^ which is also not normalizable. The surprising implication is that 
the spectrum of A is continuous. Note that the function ai is always real while /3i is 
always purely imaginary (lim^^7r/2 A = icosO). The exponent of fx is thus always real, 
implying that f\ is uniquely determined. 

4-2.2. Self-adjointness of A In the previous section we constructed a symmetric 
operator A with respect to the scalar product of Ha, i.e. A = with domain 
X'(y4) C V{A^). In this subsection we give a possible domain for A and check if there 
exists a self-adjoint extension of A. 

Definition 5. In analogy with f38l define the space CAP{R) of the (uniform) almost 
periodic functions* such that its completion is the Hilbert space L'^(IIb)- The Sobolev 
space H^(Rb) is given by the completion of the space of trigonometric polynomials 
Trig(R) in the Sobolev norm ||/||^i = ll/llL2(g^^) + II/'IIl2(Rb)' ^f^i^s) consists 
of all almost periodic functions f G CAP(R) such that f G CAP(R). 

Let the differential operator p := —i-^ on L'^{Rb) have the domain of definition 
Trig(]R). Then its closure has the domain H^{Rb)- The adjoint operator to p on L^(IRb) 
has also the domain H^{Rb) and coincides with p^ on it. Since p = p"^, p is essentially 
self-adjoint on Trig(R) [38l|39]. 

Since every almost periodic function f{x) is bounded a necessary condition for the 
inverse f~^{x) to be almost periodic is that miuj, |/(x)| ^ 0. It follows that sin~^ ^ is 
not an almost periodic function. We thus define the domain 

V{A) := {p G H\RB)®H\U{l))\p{kTX,e) = = p'^k-K^O), k G Z}, (52) 

which, according to [1I1I12], is dense. Any function ip G Vi^A) removes the pole caused 
by sin~^ ^, i.e. we require that lim^^fc^ (^(^)(sin^)~^ = and Imx^^kir ^' = 0. 
k & li. On the other hand, thanks to sin 6' in front of the differential operator id/ 89, 
the boundary term of an integration by part is automatically annihilated so that no 
boundary conditions on 6 have to be imposed. Moreover the deficiency indices n± for 
A are defined by 

n± := dimker(yl''' =]= i). 

The solutions to this equation do not lie in P(A"*") such that n±. = 0. It follows that the 
operator A is essentially self-adjoint. 

4.2.3. Quantization of p^ 3 The eigenfunctions of fP' 3 can be obtained by applying the 
same procedure on the symmetrized operator 

- Vcosd d I 

B := -ismO— -i —VcosO 
ot, sm 4 oti 

* An almost periodic function f{x) is uniformly continuous for a: g M and bounded [40] . 
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The eigenfunctions g\{i, 
M2 



Vtan(e/2) 



are given by 



cos(^/2)a2 + A 



(53) 



n(J/2)l siiiS 



fi2(cosetan({/2)), 



where 



^2 



/32 



cos^ ^ 



and 



cos^ 



cos^ 6* 

and ^2 is any C^-function that can be determined by boundary conditions. While the 
function j32 is always purely imaginary the function a2 is only real when cot^ ^/2 < cos^ 9. 
This means that g\ is not uniquely determined. We can write g\ as 

k2 



Vtan(e/2) 

with the logarithm is defined by lnF2 = LnF2 + 27rm, where n & Ij and Ln is the 
principal value of the logarithm. Inserting this solution into the eigenvalue problem 
Bgx = Xgx it can be shown that there is only a solution for n = 0. The eigenfunctions 
gx are not normalizable since the integral of 1/| tan(,^/2)| over one period is not finite. 

As for fx we are led to the conclusion that the spectrum of 13 is continuous. We can 
construct a dense subspace along the lines described in Section 14.2.21 the only difference 
being that gx has poles at ^ = 2fc7r and 6 = {2k + l)vr/2 whereas fx has poles at 
{ = (2fc + l)7r, keZ. 





0.01 



Figure 6. Absolute value of the eigenfunctions f\{(,, 0) (left panel) and g\{£,, 9) (right 
panel). The black thick line is the eigenfunction for A = 1, 6* = 1 and the black dashed 
line for A = 2, 6* = 1. 



4-3. Quantization: 2. Possibility 

In order to quantize the triad components i we replaced the configuration variables 
Qi with sm.{5iQi) / 5i. The question we may ask is to what extend this substitution 
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Figure 7. Absolute value of the eigenfunctions fx{£,, 9) (left panel) and g\{£,, 0) (right 
panel) . The black thick line is the eigenfunction for A = 1 , ^ = 2 and the black dashed 
line for A = 2, ^ = 2. 



changes the eigenfunctions. We define the symmetrized operator A2 quantized without 
the substition of Qi as 

- ^ d sin 6 d i cos 6 

The solution to the eigenvalue problem A2fx{C,, 9) = Xfx{C,, 6) is given by 

/A(e, e) = exp cos 6) v^r(log(e sin 6)), 

We see that the eigenfunctions are not almost periodic in ^. However we can choose the 
function V such that disappears, i.e. we set 

r = A/i exp ^~ 2 ^^^^^ ^^'^ ^) 
where Mi is a constant, such that the eigenfunctions to A2 are given by 

(54) 

The above eigenfunction is almost periodic in ^ but fails to be normalizable on 
L^(Rs) ® LP'{U{1)). As in the preceding section the spectrum of A2 is thus continuous. 
Note that the eigenfunction is constant in the non-diagonal limit Q ^ n j2. 
Similarly the eigenfunctions of the symmetrized operator 

- d cos Q d i sin Q 
B2 = —i sin 9— i — — I 

are given by 

^\ ,^ exp (— zA^sin^) , ^, 

gx{i, 9) = J\f2 ' -■ (55) 

vcoso' 
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The diagonal limit ^ — of (^a is just a constant function such that fP' 3, i.e. the 
expectation value (p^ 2) measures the 'diagonality' of the torus and {p^ 3) its departure. 
Once again, the eigenf unctions g\ fail to be normalizable on the Hilbert space such 
that the spectrum of B2 is continuous. Also note that contrary to the first method the 
commutator between both operators vanishes: 

[^2,^2] =0. 



4--4- Volume Operator 

The classical expression for the volume of V is given by 



V{V) 




(fx. 



Inserting the definition of the homogeneous densitized triad (fT71) we get: 
V{V)- 



|pi 1 (p2 3 -P^ 3^^2)1 (56) 

The factor t depends on the specific form of the torus and is equal to one if the torus is 
cubic such that we recover the usual expression in this limit (see e.g. Eq. (4.5) in |31j). 
Using the classical solution of the Gauss constraint we get the following expression for 
the physical volume of the torus: 



or in terms of the new variables 
V(^) 



'2 3 2 '^2V3 



h ^P^ 2 



^2^3 



1/2 



(57) 



4.4- i- Quantization of the Volume Operator according to 1. Method To perform a 
quantization of the volume operator we insert the definitions fjl3l) into ViV). Despite 
the fact that we know the eigenf unctions of the operators p^ i it is not straightforward 
to give the eigenf unctions of the volume operator V because, as explained in Section l4rT| 
they do not necessarily commute. Thus the difficult task is to determine the spectrum 
of the operator 



2J:= 



2P^ 3 



P^ 2P^ 3 



COS 



92 



sm.Q2 sinQ3 dQ"^ 
sin 6 d"^ 



— cos 



92 



+ 



dQ2dQ:, 
sin 6 (9^ 



(58) 



sin Q2, d<ddQ2 sin Q2 dQdQs ' 
However, this operator is not symmetric on Ti^ . Let us define the symmetric operator 

1 
2 



5J+ + 2J 
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- c ' If COS 6 „ sin0 d cos 6 d cos 9 d 



2 \ sin Q2 sin sin (^2 sin <96 sin dQ2 sin (^2 t^Qs / 

This operator is rather complicated and no analytic solutions to the eigenvalue problem 
could be found. 



4-4-2- Quantization of the Volume Operator according to 2. Method In this subsection 
we consider the quantization of V as described in Section 14.31 where the commutator 
between j and p"^ j vanishes. This fact simplifies dramatically the analysis because 
the (generalized) eigenvalue problem can now be written in terms of products and sums 
of the eigenfunctions of the p^ i. Let us define 

where f-y{Q2,6i), gj{Q2,6i), f!y{Q3,d2) and g!y{Q3,d2) are the (generalized) eigenfunc- 
tions of p^ 2, 3, 3 and p^ 2 respectively given in Section Furthermore we denoted 
the eigenfunctions of p^ 1 by A/ai := (Qi|Ai). Since we have 



( /a22 5'A23 ) ( /733 fl'732 , 



oc 



ex p {iQ2{\22 cos^i - A23 sin 9i)) exp (2(^3(733 cos6'2 - 732 sin 612)) 
Vsin 9i cos6'i V sin 62 cos 62 



we see that 7ai,A22,A23,733,732 normalizable in . The generalized eigenvalue 

problem is thus given by 

VT\j^,A22,A23,733,732 ~ "^Ai , A22 , A23 ,733 1732 ["^''/-'] 

= 7^^^^Pl\/fc|Al(A22733 - A23732)|Tai,A22,A23,733,732 M (59) 

for if G V{V). 



4-5. Quantum Gauss Constraint 

In Section 13.4.21 we computed the classical Gauss constraint for a Bianchi type I model. 
In the open case the elementary variables can always be diagonalized such that both the 
diffeomorphism and Gauss constraints are automatically satisfied. In the closed model 
this is not the case anymore so that a quantization of the constraints is mandatory. 
Since in Bianchi type I models the diffeomorphism constraint is proportional to the 
Gauss constraint we only need to quantize and solve the latter. However, contrary to 
the diffeomorphism constraint the Gauss constraint can be quantized infinitesimally. 
A gauge transformation of an SM(2)-connection is given by 

A' = X~^AX + X-^dX 

where A : S i-^ SU{2). Infinitesimally we can write this equation as 

Aj. ^ A': = A: + dae' + ,,eM^ + 0{e^). 
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The classical Gauss constraint ensuring 5'[/(2)-invariance is given by 
G(A) = - / £xE''DaA^ 

where DaA^ = daA^ + e-' ki^a-^'' is the covariant derivative of the smearing field A-' . The 
infinitesimal quantization of this expression yields an operator containing a sum of right 
and left invariant vector fields over all vertices and edges of a given graph a. This 
operator is essentially self-adjoint and can, by Stone's theorem, be exponentiated to a 
unitary operator t/^ defining a strongly continuous one-parameter group in 0. Usually, 
in order to find the kernel of the Gauss constraint operator one restrict the scalar 
product on Tiaux to the gauge-invariant scalar product on Ti.^-^. This Hilbert space is 
a true subspace of Tiaux since zero is in the discrete part of the spectrum of the Gauss 
constraint operator. 

We saw in Section 13.51 that thanks to the symmetry reduction two of the Gauss 
constraints are automatically satisfied. While the nonvanishing Gauss constraint (l30l) 
is still a complicated function in and p"^ j it simplifies to Equation fl35l) after the 
canonical transformation. A quantization of this expression is then given by 

Gi = Pq^ — Pg^ ■ 

Since the eigenstates of the momentum operators Pg^ are the strict periodic functions 
satisfying Equation (HOj) the action of the Gauss constraint on k) is given by 

which vanishes if 

h = k2. 

We can thus introduce a new variable := + ^2 such that the algebra As given by 
Equation fl38|) reduces to the invariant algebra A^g"^ generated by the functions 

X exp QzAiQi + ^iX2Q2 + ^^AsQa + ikO^ . (60) 

A Cauchy completion leads to the invariant Hilbert space Hf^^ = Tif^ x Ti^i. A 
comparison with Ti^ shows that we 'lost' one Hilbert space 7^51 by solving the quantum 
Gauss constraint. Furthermore, instead of two momentum operators conjugated to 61 
and 02 we have just one momentum operator conjugated to defined by 

The eigenstates of all momentum operators are given by 

l/i, k) := \fii, fi2, f^3, k), 



where /c G Z defines the representation of U{1). 
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5. Conclusions and Outlook 

In this paper we studied how a torus universe affects the results of LQC. To do so 
we first introduced the most general tori using Thurston's theorem and found that six 
Teichmiiller parameters are needed. We construted a metric describing a fiat space but 
respecting the periodicity of the covering group used to construct the torus and used it to 
derive a gravitational Hamiltonian. We studied the dynamics of a torus universe driven 
by a homogeneous scalar field by numerically solving the full Hamiltonian and saw that 
its form only remains cubic if all off-diagonal terms vanish. The Ashtekar connection 
and the densitized triad for a torus were then derived for both the most general and a 
slighty simplified torus. The reason for this simplification was that a simple solution 
to the Gauss constraint could be given. We also derived the Hamiltonian constraint in 
these new variables and showed that it reduces to the standard constraint of isotropic 
LQC in case of a cubical torus. 

The passage to the quantum theory required a canonical transformation so as to 
be able to write the holonomies as a product of strictly and almost periodic functions. 
A Cauchy completion then led to a Hilbert space given by square integrable functions 
over both and U{1). However the drawback of the canonical transformation is a 
much more complicated expression for the components of the densitized triad containing 
both the momentum and the configuration variables. Following the standard procedure 
of LQC we substituted these configuration variables with the sine thereof and were 
able to solve the eigenvalue problem analytically. Surprisingly it turned out that the 
(generalized) eigenfunctions of the triad operators do not lie in the Hilbert space, i.e. 
the spectrum is continuous. On the other hand we were also able to find almost 
periodic solutions to the eigenvalue problem of the triad operators without performing 
the substitution just described, but once again these eigenfunctions do not lie in the 
Hilbert space. The reason why both ways lead to a continuous spectrum is the non- 
cubical form of the torus, for if we set the angles 6*1 2 = in Equation fl32l) the triads 
correspond to the ones obtained in isotropic models. Furthermore we were able to find 
the spectrum of the volume operator for the second case because, contrary to the first 
case, it is a product of commutating triad operators. 

Although we gave a couple of numerical solutions to the classical Hamiltonian we 
didn't consider its quantization. The constraint describing quantum dynamics of a torus 
is given by inserting the holonomies ( 137|) into Thiemann's expression for the quantum 
Hamiltonian [13] 

Contrary to LQC and LQC we saw that the spectrum of the volume operator of a torus 
is continuous. It would thus be very interesting to know how far Cgrav departs from the 
usual difference operator of LQC. Furthermore, whether a quantization of the torus a la 
LQC removes the Big Band singularity needs also to be addressed, especially since we 
saw that many characteristics of both LQC and isotropic LQC are not present in this 
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particular topology. 

In this work we only considered the simplest closed fiat topology but there are many 
other closed topologies. As we saw there are eight geometries admitting at least one 
compact quotient. For example there are six different compact quotients with covering 
E^, namely T^, T^/Z2, T^/Za, T^/Z4, T^/Zg and a space where all generators are screw 
motions with rotation angle 7r/2. It would be interesting to know how these discrete 
groups Z2,3,4,6 affect the results of this work, especially since the last five spaces are 
inhomogeneous (observer dependent) [ii] . 
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Appendix A. Fundamental Domain of the Torus 

In two dimensions the upper half-plane H is the set of complex numbers H = {x+iy \ y > 
0; X, ?/ G R,}. When endowed with the Poincare metric 

ds"^ = {dx^ + dy^^jy^ this half-plane is called the Poincare upper half-plane and is a 
two-dimensional hyperbolic geometry. The special linear group SL{2, R) acts on H 
by linear fractional transformations z i— > {az + h)/{cz + d), a,b,c,d e IR, and is an 
isometry group of H since it leaves the Poincare metric invariant . The modular group 
S'L(2, Z) C S'L(2,R) defines a fundamental domain by means of the quotient space 
H/SL(2, Z). This fundamental domain parametrizes inequivalent families of 2-tori and 
can thus be identified as the configuration space of the two-dimensional torus. Since 
we consider a three-dimensional torus with six independent Teichmiiller parameters (see 
Equation ([2])) we need a generalization of the Poincare upper half-plane [l5l |16] to a 
six-dimensional upper half-space. 

Definition 6. A fundamental domain D for SL{3, Z) is a subset of the space 
^3 := {A E Mn{3, R) | A'^ = A, A positive definite} which is described by the quotient 
space ^s/SL{3, Z). In other words, if both A G ^3 and A[g] := g^Ag, g G SL{3, Z), 
are in the fundamental domain then either A and A[g] are on the boundary of the 
fundamental domain or g = id. 

Since ^3 is a subspace of the six-dimensional Euclidean space (there are six 
independent matrix elements for A G ^3), the generalization of the Poincare upper 
half-plane is now a six-dimensional upper half-space := {(ai, . . . , ag) G | ag > 0} 
upon which the group SL{3, R) acts transitively. To identify with an upper half- 
space we introduce the Iwasawa coordinates such that VA G .^3 there is the unique 
decomposition: 
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with Xj, Uj G R- with Y\yi = 1. The geometry of the upper half-space is given by the 
GL{3, R)-invariant hne element 

Note that the Ricci scalar of the metric (lA.ip is constant and negative. 



In order to give a parametrization of the fundamental region we use Minkowski's 
reduction theory [17], which tells us that for a metric hi^j the following inequalities must 
be satisfied: 



hi^i < /ij+i^j+i, i — 1,2,3 
1 
2 

Since the metric (jlj) is invariant under the map 03 ^ — » —03 ^ we can define the upper half- 
space as = {(oi ^, 02 ^, (32 ^5 ^5 ^5 % ^) € I 03 > 0}, where we have identified 
the element with ^. In our parametrization (jll) we therefore obtain the fundamental 
domain D delimited by the inequalities: 

1 1 1 

^2 < -ai 

1 1 1 
03 < 2*^1 

a2 '03 ' + 02 '03 ' < ^ ((^2 + (a2 

The first inequality tells us that the length of the generators of the torus are ordered: 
llc^ill < ||a2|| < ||fl3||- However, starting with such an ordered triplet does not necessarily 
imply that the order is preserved by dynamics. Thus we think that it may be more 
appropriate to choose the equivalent representation of the configuration space given by 
C = R®. Otherwise, we would have to rotate the coordinate system every time the torus 
leaves the fundamental domain. Note that similar results have also been obtained in 
M-theory, where one considers a compactification of the extra dimensions on T" (see e.g. 
[miin]). However, the situation is different in string theory where one really integrates 
only over the fundamental domain, e.g. Z(T") = J^drZi^r). 

Appendix B. The Torus Universe in Iwasawa Coordinates 

In this appendix we use a parametrization of the torus using the Iwasawa coordinates 
which are more apt to describe the asymptotic behavior of the metric [50]. It is important 
to understand the role of the off-diagonal terms in the metric (jll) and to know what 
happens near the singularity and at late times. The metric can be decomposed as 
follows: 

h = X^V^Af, (B.l) 
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where 
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An easy calculation shows that Equation (IB.ip can be transformed into Equation (jlj) 
with Til = 02 ^/ai^, n2 = a^^/ai^, = 03^/02^, Zi = — Inoj* (no summation) jj. The 
analogue to Equation (j3]) is now given by 

3 



(B.2) 



i=l 



The Iwasawa decomposition can also be viewed as the Gram-Schmidt orthogonalization 
of the forms dx"". 

3 



i=l 

where the coframes 6^ are given by 

=Arjdx\ 

Analogously, the frames Cj dual to the coframes 6^ are given by the inverse of A/l ' 



' dx"- 



Since the determinant of the matrix J\f is equal to one the basis given by the coframe is 
orthonormal. Note that this is different from the construction in Section 12.21 

In order to determine the asymptotic behavior of the off-diagonal terms we follow 
the analysis in [50]. The metric h being symmetric, we automatically know that its 
eigenvalues are real. We call these eigenvalues 1 < « < 3 and ai < a2 < a^, 
in analogy to the diagonal Kasner solution (see Section 13.4.11) and construct a metric 
hxit) = diag(t^"') by means of a constant matrix L 

\ 



h{t) = L^hK{t)L, L 



h h ^3 

mi 7712 

ri 72 73 j 

With these relations we can obtain the evolution of the Iwasawa variables. For example 
we have 

_ ^^"'^1^2 + t^"'mim2 + t2°3^ir2 

In [50] it was shown that the asymptotic behavior t 
the Iwasawa variables is given by 

k h hms - h7ni 



rii 



n2 



k 



lim2 - hrrii ' 



0"^ of the off-diagonal terms of 
(t-0+). 



jJ For simplicity we assume that all diagonal scale factors ' are strictly positive. However the 
nondiagonal scale factors can be negative or zero. 
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which means that the off-diagonal terms freeze in as we approach the singularity. We 
can also calculate the other limit t ^ oo and obtain e.g. 

rii -> — , 722^—, ^3 -> , [t^oo). 

ri ri mir2 — m2ri 

We have checked this result numerically, which can be seen in Figure [2] where the gray 

line (a2 ^) converges for t —> cxd toward the solid line (ai i.e. rii — > const. However, 

the limit t — 0"*" could not be checked due to the numerical instability of the solutions 

when approaching the singularity. 
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